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We propose a model that can account for the experimentally observed phase behavior of DNA- 
nanoparticle assemblies (R. Jin et al., JACS 125, 1643 (2003); T. A. Taton et al., Science 289, 
1757 (2000)). The binding of DNA-coated nano-particles by dissolved DNA linker can be described 
by exploiting an analogy with quantum particles obeying fractional statistics. In accordance with 
experimental findings, we predict that the phase-separation temperature of the nano-coUoids in- 
creases with the DNA coverage of the colloidal surface. Upon the addition of salt, the demixing 
temperature increases logarithmically with the salt concentration. Our analysis suggests an exper- 
imental strategy to map microscopic DNA sequences onto the macroscopic phase behavior of the 
DNA-nanoparticle solutions. Such an approach should enhance the efficiency of methods to detect 
(single) mutations in specific DNA sequences. 

PACS numbers: 



Colloidal particles that can be selectively linked by 
specific sequences single-stranded DNA, represent an en- 
tirely novel class of complex liquids. With these par- 
ticles, it becomes feasible to "design" multicomponent 
mixtures, where the attractive interaction between every 
pair of species can be controlled independently. 

During the past few years, the unique molecular recog- 
nition properties of DNA have been exploited to cre- 
ate self- assembled nano-structures and nano-devices with 
a remarkable degree of control (see e.g. Ref. P]). 
Self-assembled nanostructures based on DNA have also 
been utilized to detect specific DNA sequences 0, 
0, 0, In nanoparticle-based DNA detection sys- 

tems, single-stranded DNA (ssDNA) "probe" molecules 
are chemisorbed onto the surface of gold nanoparticles 
(Fig. These probe ssDNA's can specifically bind to 
a dissolved "target" ssDNA and thus detect its presence 
with high sensitivity and selectivity: at high enough con- 
centration, the target ssDNA strands induces a sharp 
demixing transition that leads to aggregation of the 
nanoparticles coated with the complementary "probe" 
strands. These aggregates can be easily detected by op- 
tical means. 

The sensitivity of the nanoparticle-based DNA detec- 
tion method is some two orders of magnitude higher than 
that of the corresponding fluorophore-based DNA-array 
scheme (sf. But, in addition, the method can be used to 
distinguish different DNA sequences that differ from each 
other by only a single base. The reason is that a single 
mismatch in the DNA sequence results in a significant 
change in the DNA-coUoid demixing temperature. The 
fiuorophore-based systems lack this level of selectivity . 

In a recent paper, Jin et al. Q have analyzed the de- 
pendence of the dissolution temperature (see note 6|) of 
DNA-linkcd nanoparticles on the DNA coverage density, 
the particle size, and the salt concentration. In the ex- 
periments 0,3 short DNA molecules (~ 20-base oligonu- 
cleotides) and small gold colloids (~ 10 — 50 nm) were 
used. The principal experimental findings are the follow- 



ing: (i) The higher the ssDNA coverage density on the 
surface of colloids, the sharper the dissolution profiles, 
and the higher the dissolution temperature, (ii) The ad- 
dition of salt stabilizes the aggregates. The dissolution 
temperature increases logarithmically with salt concen- 
tration. The higher the salt concentration, the larger the 
total fraction of DNA-nanoparticle aggregates. 

Below, we propose an explanation for these experi- 
mental findings. We argue that the origin of the sharp 
phase-transition profiles is the entropic cooperativity of 
the DNA-nanoparticle network. Upon cooling, the sys- 
tem undergoes liquid-liquid phase separation. The dense 
"liquid" phase is strongly cross-linked and behaves as a 
solid gel. 

We first consider a simple model for a binary mixture 
of DNA-coated colloids. Each colloidal species is covered 
with specific ssDNA molecules, in such a way that all 
possible pairs of colloids can be bound by three types of 
complementary ssDNA linkers: viz., AA, BB, and AB. 

For the sake of simplicity, we use a lattice model to an- 
alyze the phase behavior of this system. Colloids of type 
A and B can occupy the nodes of a three-dimensional 
lattice with coordination number q. The system is grand 
canonical with respect to all the components: /ij^ and 
fig are the chemical potentials of colloids of species A 
and B, respectively; and (J-aa, I^bb, and fiAB a-i'e the 
corresponding chemical potentials of linkers. First, we 
consider the situation where only linkers of type AB are 
present (this corresponds to the experimental conditions 
of Refs. 0, For this system, the grand canonical 

partition function is: 
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where Na and Nb are the number of colloids A and B, 
respectively, n is the number of nearest neighbor AB col- 
loidal pairs for a given realization of the grand canonical 
ensemble, and g(N a-, Nb, n) is the total number of possi- 



ble configurations with n nearest neighbor AB pairs for 
a given values of Na and Nb ■ The sum with respect to n 
extends over all values of n consistent with the fact that 
there are Na and Nb colloids A and B present. The 
number of links between a given pair of colloids depends 
on the surface density of probe strands on the A and B 
colloids. If there is only one pair of probe strands per con- 
tact, then this pair can accommodate at most one linker 
and the bound linkers obey Fermi-Dirac statistics. In 
contrast, when there are (infinitely) many probe strands 
per contact, then the binding of linkers is determined 
by Bose-Einstein statistics 0. In general, the maximum 
number of linkers per contact is a number M > 0. For 
cxD > Af > 1, the bound linkers obey fractional statistics. 
In practice, the local coverage of probe molecules fluctu- 
ates . We assume that the maximum number of linkers 
per bond, obeys Poisson statistics. Moreover, we assume 
that the values of M for different colloid pairs, are uncor- 
related. If we average over all possible values of Af, we 
obtain the following expression for the grand partition 
function of linkers (for a given AB colloidal pair): 
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where M is the average value of Af, and z = e t . 
Eq/j (T) is the temperature-dependent binding free energy 
of a double-stranded DNA (dsDNA) molecule connecting 
a given pair of aP colloidal particles. By tuning Af, we 
change the effective fractional statistics of the linkers. 
The Bose-Einstein limit {Qab = 1/(1 — ^)) is recovered 
when Af ^ 1. The Fermi-Dirac limit is only recovered 
when Af is fixed, and equal to one. The physics of DNA 
melting and its effect on the phase behavior of the DNA- 
colloid system is described exclusively by Eq/j. 

It is straightforward to generalize the model (Eq. Q), 
to the case where all three possible types of linkers are 
present. This results in multiplicative factors in Eq. 
that account for all possible configurations of all link- 
ers: QYb" Qaa Qslf^ where Qaa and Qbb are defined 
analogously to Qab, Eq. |(2Jl. 

The expression given by Eqs. can be interpreted 

as the partition function of a three-state spin model with 
Hamiltonian 
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Ui = (1, 0), (0, 1), or (0, 0) if A, B, or a vacancy (solvent) 
occupies site i; and fj = {fJ-A,tJ-B)- The square, symmet- 
ric interaction matrix is given by J"^ = —Tin Qa0- We 
note that the effective attraction between the colloids in- 
duced by linkers is dominated by the entropy associated 
with the number of different ways of distributing L link- 
ers over M bonds: in a dense phase, there are simply 
more bonds. Similar, entropic mechanisms are responsi- 
ble for phase separation in binary hard-core mixtures Ig , 
polymer- microemulsions 9], and microemulsions |lCllll| . 
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FIG. 1: Computed temperature dependence of tlie fraction 
of the ssDNA-coated colloids in tlie fluid (dilute) phase of a 
binary colloidal mixture (species A and B) and added ssDNA 
linkers. Above the demixing temperature, the volume frac- 
tions of both the A and B colloids are equal to 20%. The tar- 
get ssDNA can form links between A and B colloids, but not 
between A and A, or B and B. At low temperatures, virtually 
all colloids end up in the dense "gel" phase. The two curves il- 
lustrate the effect of a mutation in the target DNA. The solid 
and dashed curves were obtained with fi = hab = 3.5, s = 10, 
g = 6, M — 5. All temperatures are relative to the melt- 
ing temperature of the perfectly complementary linker (solid 
curve). The dsDNA melting temperature for the dashed curve 
(corresponding to a base-pair mismatch) is 10 % lower than 
Tm for the perfectly matched linker. The arrow indicates the 
optimal stringency temperature. 



The Hamiltonian for a three-component lattice gas 
model (Eq. ©) has been studied extensively 0|. The 
mean- field expression for the free energy is: 



/ 



T4)A In 

r(i - 



4>A 



f T(j)B In 4>B 
-0s)ln(l- 

^^^B + 



qj 



AB 



(4) 



where (pA and (j)B are the average volume fractions of 
A and B colloids. The equilibrium phase behavior of 
the system in the mean-field approximation follows from 
the analysis of this free energy 0|: depending on the 
strength of the interaction between the colloidal species, 
the system can be in a homogeneous state, or separate 
into two, or even three, coexisting phases. It is straight- 
forward to generalize the model to any number of col- 
loidal species and corresponding linkers. 

We first apply the above model to analyze the exper- 
iments of Refs. I^Q. In those experiments, only AB 
linkers were present in solution, and hence J^^ = J^^ = 
0. The free-energy difference between the helix and coil 
states of a dsDNA molecule, eAB{T) = e(T) constitutes 
an input to the interaction potential. We adopt the sim- 
ple form e(T) = s {T — T^), where is the melting 
temperature of a dsDNA molecule and s « 10 ,14]. The 
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dissolution curves of the DNA-linked colloidal aggregates 
are shown in Fig. ^ The fraction of colloids in solution 
is plotted as a function of temperature. Above a criti- 
cal temperature Tc, all colloids are in solution. Consis- 
tent with the measurements of Ref. Q, the dissolution 
temperature of the aggregates increases with increasing 
surface coverage of ssDNA molecules (represented by the 
linker occupation number M) on the colloids. In the 
"Bose-Einstein" limit (M ^ 1), Tc saturates at a value 
= in[i-cxpt-8/g)]+. ' as follows from Eq. Q in the 
case of equal initial volume fraction of colloids A and B 

The fact that the dissolution of the DNA-linked col- 
loidal aggregates occurs in a narrow temperature range, 
makes it possible to distinguish DNA strands that dif- 
fer in only a single base-pair ,3J. This is achieved by 
selecting a temperature ("optimum stringency temper- 
ature" ) where a perfectly matched target results in the 
aggregation of a high proportion of the nanoparticles, 
while the same amount of ssDNA with a single mismatch 
causes only a small fraction of the nano-particles to ag- 
gregate. As can be seen from Fig. ^ our model repro- 
duces the observed behavior: For the model calculations 
shown in Fig. ^ the optimum stringency temperature 
(indicated by the arrow) is the dissolution temperature 
(T on the binodal) for the mismatched-DNA-linked col- 
loids, T = 1.273 (dashed curve) (in the units of T,„ for 
the perfectly complementary dsDNA), and we assumed 
that T„i for the mismatched DNA is 10% lower com- 
pared with Tjn of the perfectly complementary DNA. At 
this temperature 100% of the mismatched and only 0.1% 
of the perfectly complementary DNA-linked colloids are 
dissolved. The selectivity is thus extremely high. 

Next, we consider the effect of the varying of salt con- 
centration on the melting of the DNA-linked colloids. 
Added salt has two effects: It changes the DNA melt- 
ing temperature and it modifies the electrostatic colloid- 
colloid and DNA-coUoid interactions. For the salt con- 
centrations, used in the experiments of Ref. JH, all 
electrostatic interactions are strongly screened [Tjl and 
hardly depend on the salt concentration. However, the 
salt concentration does affect the melting temperature of 
the DNA. To estimate this effect, we make use of the 
known properties of poly-electrolytes, in particular, the 
theoretical prediction [l3| that the melting temperature 
of the DNA helix-coil transition is proportional to the 
logarithm of the salt concentration. This result has been 
verified experimentally (see e.g., Ref. 0|). A more de- 
tailed analysis of this result in the context of the present 
work will be presented in Ref. [l^ . Here we simply 
use the fact, that for experimentally used salt concen- 
trations, the DNA melting temperature varies with salt 
concentration as 0, 0|: Tm = [1 -t- alogc], where 
is the DNA melting temperature at a reference state 
(with a salt concentration Iq = 0.1 M), a « 0.05, and 
c = I/Iq. We then obtain the following expression for z: 
z — exp[ ^^^ j° ^"^"^^ — s], where /i and T are expressed 
in the units of T'^-^. This expression reproduces the ex- 
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FIG. 2: Selectivity of the phase behavior of a ternary sus- 
pension consisting of j4, B, and C colloidal species and one 
type of ssDNA linker that is perfectly complementary to AB 
pair and slightly mismatched with respect to AC pair (see 
text). The solid (dashed) curves show that the presence of 
AB linkers induces a large relative concentration differences 
(0fl — 4'c)/{4'B + (j^c) of B and C colloidal species in the 
dense (solid curve) and dilute —{(j>B — <l>c)/{4>B + 4'c) (dashed 
curve) phases, respectively. Above the demixing tempera- 
ture, the concentration of all colloidal species is the same: 
'j>A = 4>B = 4>c ~ 0.17. The parameters used in to compute 
this figure are: /i = 2.5, s — 10, q = 6, M — 3. T and fi 
are expressed in the units of Tm of the linker that binds AB. 
The melting temperature of the linker connecting AC was as- 
sumed to be only 5% lower than T,n- The arrow indicates the 
optimal stringency temperature. 



The model described above suggests a novel approach 
to DNA screening (mutation analysis). As mentioned 
above, colloids functionalized with oligonucleotides have 
been used as probes to differentiate perfectly complemen- 
tary targets from those with single-base mismatches (see 
Ref. Q)- The drawback of the existing method is that 
the thermodynamic properties of DNA-linked nanopar- 
ticle assemblies are sensitive to the parameters of the 
system, such as the grafting DNA density, the concen- 
trations of DNA linkers and colloids, the monodispersity 
of colloids, and the pH and ionic strength of the solvent. 
As it is essential to maintain precisely the same experi- 
mental conditions when performing the measurements on 
different target DNA strands, mutation-screening exper- 
iments are rather demanding. 

Here we propose a method that should be relatively 
insensitive to unavoidable variations in the experimental 
conditions. To illustrate our approach, consider a sys- 
tem consisting of three colloidal species A, B, and C 
(see Fig. 121. Species B and C are grafted with ssDNA 
molecules that have similar sequences and only different 
by a mismatch (either a single-base or multiple-base) . We 
assume that the target ssDNA is either perfectly comple- 
mentary to AB sequence, or mutated, and thus perfectly 
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complementary to AC sequence (the method is easily 
generalized to a larger possible number of mutations). 

The addition of the target ssDNA, induces a two-phase 
separation in the system. Importantly, the composition 
of the dense aggregate (and, by implication, the com- 
position of the remaining solution) depends on the na- 
ture of the target ssDNA. The characteristic mean-field 
phase diagram for this case is shown in Fig. |21 in terms 
of the concentration differences (selectivity plot) between 
B and C colloids. The target is perfectly complementary 
to AB. The melting temperature of the target DNA is 
assumed to be some 5% higher for a perfectly matched 
sequence than for a mismatched sequence. The initial 
concentrations of A, B, and C colloids are chosen to be 
equal. Our main prediction is that the sequence of the 
target can be determined simply by monitoring the con- 
centration differences of the species B and C in either a 
dense (solid curve) or dilute (dashed curve) phase. In the 
case shown in Fig. |21 at the optimal stringency tempera- 
ture T « 1.165 (in the units of of the AB linker), the 
difference between the concentrations of B and C colloids 



in the aggregated phase is as large as 30%. If the colloids 
fluorescently labeled, even small differences in concentra- 
tions can easily be detected. In this way, a single experi- 
ment allows us to distinguish between the perfectly com- 
plementary and mismatched target. We stress that the 
selectivity of the DNA screening can be tuned by varying 
the surface coverage density M. The method should be 
robust to moderate variations in the experimental condi- 
tions. It is straightforward to generalize the technique to 
discriminate between a larger number of possible muta- 
tions, simply by increasing the number of colloidal probe 
species. 

In essence, the method allows us to map the micro- 
scopic base-pair sequence of the target ssDNA onto the 
macroscopic phase behavior of the DNA-linked nanopar- 
ticle solution. 
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